Quantum magnetism with multicomponent dipolar molecules in an optical lattice 
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We consider bosonic dipolar molecules in an optical lattice prepared in a mixture of different 
rotational states. The 1/r 3 interaction between molecules for this system is produced by exchanging 
a quantum of angular momentum between two molecules. We show that the Mott states of such 
systems have a large variety of quantum phases characterized by dipolar orderings including a state 
with ordering wave vector that can be changed by tilting the lattice. As the Mott insulating phase 
is melted, we also describe several exotic superfluid phases that will occur. 
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In ultracold physics, systems with long-range dipolar 
interactions have recently attracted considerable atten- 
tion both theoretically and experimentally (for a recent 
review of ultracold dipolar molecules see |l| and refer- 
ences therein). For atoms, dipolar interactions come from 
their magnetic moments and become important for large 
electronic spin Q ■ Recent experiments demonstrated the 
relevance of such dipolar interactions for the expansion of 
Cr atoms from the BEC state Q. On the other hand, for 
heteronuclear molecules, dipolar interactions arise from 
their electric dipole moments. Recent experiments have 
succeeded in trappingand cooling several types of het- 
eronuclear molecules [ll Li la, l1 • I n a state with a well- 
defined angular momentum, molecules do not have a 
dipole moment. However, when an external electric field 
is used to polarize the molecules, dipolar moments can be 
induced. There has been considerable theoretical effort 
to study the resulting dipole interactions and many - body 
physics associated with such systems [a 111 111! LLU, 1 121 ? ] . 

In this Letter, we consider an alternative mechanism 
for obtaining the 1/r 3 dipolar interactions, and the im- 
portant concomitant directional character. Namely, we 
investigate a mixture of heteronuclear dipolar molecules 
in the lowest (Af = 0) and the first excited (Af = 1) rota- 
tional states. For such a system, the origin of the long- 
range interaction is the exchange of angular momentum 
quanta between molecules. We demonstrate that when 
loaded into an optical lattice, such mixtures can real- 
ize various kinds of non-trivial effective dipolar spin sys- 
tems with anisotropic, long-range interactions. Several 
approaches for realizing spin systems using cold atoms 
have been discussed before, includ ing bosonic mixtures 
in optical lattices in the Mott state Ha . FLU Il5l H(| , inter- 
acting fermions in special lattices |l7j . and trapped ions 
interacting with lasers [l^. The system we consider has 
the practical advantages of the high energy scale for spin- 
dependent phenomena (set by dipolar interactions) and 
the new physics associated with the long-ranged nature 
of the dipole interactions. Experimental realization of 
the system will give insight into several open questions in 
condensed matter physics including competition between 
ferro and antiferroelectric orders in crystals 0, and 
systems with frustrated spin interactions |2lj . 



Consider the system that contains bosonic molecules 
in the lowest (AC = 0,AC 2 = 0) and first excited (AC = 
1, Af z = —1,0, 1) rotational states where we let and 
*_i 1 create these respective states. We will often use 
the change of basis 4 = (t\ + *Li)/V2, *J = -i(t{ - 
£_ 1 )/>/2j and t\, = t\. To describe molecules in an optical 
lattice we use the one-band Hubbard type effective model 



W — Wkin + ?^Hub + ^dip- 
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The first term on the right hand side of JTJ is the 
kinetic energy from nearest-neighbor hopping 7ikm = 



- J E(y> (4 s i + ****** + h - c -J ■ Operators s] and t] a 
create molecules on site i (here and after the summa- 
tion over repeating indices a = x, y, z is implied). The 
last term in \$ describes the dipolar interaction between 
molecules from different sites 
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where R^ are lattice vectors, Bij a is the a-component of 
the unit vector along R^ — Rj, and parameter 7 equals 
2cP/3, where d is the value of the dipole moment associ- 
ated with the Af = — > Af = 1 transition, and di is the 
dipole moment operator at site i. The a-component of 
the operator is written as 



t-j -2iB e t 1 A HB r t 
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where we absorbed the energy difference between the ro- 
tational levels Eflf = i — £V=o = 2-B e into the time depen- 
dence of the t operators. Since the rotational constant 
B e is considerably larger than any other energy scale in 
the system, we assume that the terms in © that oscil- 
late at frequencies ±4-B e average to zero. This forces the 
number of molecules in the Af = and Af = 1 states to 
be independently conserved. Then (J2J) reduces to 



^dip = 2 E 
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The second term on the right hand side in is the 
Hubbard on-site interaction. For two s molecules in 
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the absence of an external electric field, the long-range 
part of their interaction potential is dominated by the 
van der Waals tail C e /R 6 originating from second order 
terms in the dipole-dipole interaction operator j^. For 
polar molecules with large static rotational polarizabil- 
ities one can estimate C§ w — d 4 /6B e . For the RbCs 
molecule (d = 0.5 a.u., B e = 7.7 x 10~ 8 a.u.) we have 
C*6 ~ 1.5 x 10 5 a.u. For molecules with smaller dipole 
moments and larger rotational constants like, for exam- 
ple, CO (d = 0.043, B e = 9.0 x 10 -6 ), the van der Waals 
interaction is comparable in magnitude to interatomic 
forces. In any case the range of the potential, which 
scales as R e = (mCg) 1 ^ 4 , is not much different from typ- 
ical ranges of interatomic potentials (for RbCs R e ss 400 
a.u.). First order terms in the dipole-dipole operator are 
also absent for two molecules with TV = 1. In this case, 
apart from a weak quadrupole-quadrupole contribution 
proportional to R~ 5 , the long-range part of the inter- 
molecular potential is given by the van der Waals in- 
teraction with a comparable Cq coefficient. Thus, the 
interactions between molecules with the same AT are all 
short ranged and, in an ultracold system, can be mod- 
eled by contact potentials. Then, averaging them over 
the Gaussian on-site wave functions gives the Hubbard 
on-site interaction. 

The interaction between s and t a molecules (without 
loss of generality we consider a = z here) is similar 
to the resonant interaction of an electronically excited 
atom and a ground state atom. For even partial waves 
the intermolecular potential is asymptotically given by 
W Z (R) = 7(1 - 3cos 2 9 Z )/2R 3 , where 9 Z is the angle 
between R and the z-axis. We consider the weakly inter- 
acting regime where the characteristic energy scale of this 
interaction, AE ~ jIq , is smaller than the Bloch band 
separation. Here lo is the oscillator length of the on-site 
harmonic confinement. Then, the two-body problem in 
a harmonic potential can be solved in the mean-field ap- 
proximation by using the pseudopotential approach (see 
|23| and references therein). Due to the anisotropy of 
W^(R) the corresponding on-site interaction energy, V z , 
can be tuned at will by changing the aspect ratio of the 
on-site confinement [23j. 

We arrive at the following expression for H.n u b'- 
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It is easy to see that (JSJ holds for arbitrary filling factors 
as long as the on-site density profiles remain Gaussian. 
However, for the purpose of this paper it is sufficient to 
consider on average one or two molecules per site. The 
ten coupling constants in Eq. 10 are difficult to control 
independently. However, in many cases not all of them 
arc relevant. For a mixture of s and t, molecules the 



relevant coupling constants are U, U z , and V z . These 
arc tunable through trap aspect ratio and/or Fcshbach 
resonance. When considering this system we will take 
U = U z and V = V z . Another example is the Mott 
insulating state with one molecule per site (all kinds of 
molecules are allowed) when the interaction energies are 
much larger than J. Then, the particular values of the 
on-site coupling constants are not important (for me- 
chanical stability it is sufficient that they are positive) 
and the state of the system is found by minimizing the 
intersite dipolar interactions. 

We now discuss the resulting phase diagram, first fo- 
cusing our attention on the Mott insulating state with 
one molecule per site. We take the variational wave func- 
tion 

i* M i) = n ( c ° s ( e )4 + M^icti) 10) (6) 

i 

where 9 describes the fraction of the molecules excited 
into Af = 1 states and ipi a is a normalized complex vector 
- 0i*a'0ia=l- Here, ipi a is the variational parameter which 
descibes the direction the dipole moment points on site i. 
This allows us to construct variational states that bene- 
fit maximally from dipolar interactions. In all cases dis- 
cussed below we verified the absence of phase separation 
by checking the eigenvalues of the compressibility matrix 
for s and t bosons [24|]. Taking the expectation value of 
the dipole operator with our variational wave func- 
tion (J5J we obtain (di a ) = sin(29)\ip a \cos(<fi a — 2B e t) 
where we have written ipi a = lipale^^ . Upon taking the 
expectation value of dipole Hamiltonian we find for 
the dipolar energy 
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When minimizing the energy in (JJJ it is important to 
keep track of the conservation laws that may be present 
for certain experimental geometries and on the initial 
preparation of the system. We will now consider sev- 
eral examples of ordering in the Mott insulating state. 
Although the dipole interaction in all cases is described 
by (7J we will see that different preparation leads to very 
different types of order. Though all discussion in this 
work will be restricted to 2d, we emphasize that there 
are nontrivial results in the Mott insulating phase for the 
Id and 3d cases as well. As the first example, we con- 
sider the square lattice in the xy-plane defined by vectors 
ai = x and &2 = y- Due to cross-terms such as sH\.st y in 
the dipolar hamiltonian we see that t x molecules can 
be converted to t y molecules and vice-versa. Thus, Nt x 
and N ty are not conserved quantities, and, consequently, 
the only conserved quantities arc N s and N tz . Now 
consider preparing this system in a mixture of Af = 
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and Af = l,Af z = 1 states. Then after the system re- 
laxes, taking the constraints into account, we must have 
fixed (N s ) = Ncos 2 (9), (N t J + (N ty ) = N sin 2 (9) , and 
(N t _) = 0. This gives the constraints on the variational 
wave function if,> iz = and \ip x \ 2 + IV'yl 2 = 1- We see 
that the dipoles are allowed to rotate freely in the xy- 
plane. For this case, the dipoles will choose to point 
head-to-tail in the direction of one of the bonds, while 
alternating in the other. Thus, it is straightforward to 
see that this gives the ordering wave vector q = (0,7r,0) 
with ipi x = e*( q ' Ri+v °) and ipiy = tyiz — where (fo is 
an arbitrary phase corresponding to a change of phase of 
the time dependent oscillations of the dipolar moment. 

We point out that this configuration is degenerate 
to the one with dipoles pointing head-to-tail in the y- 
direction. 

As the next example in two dimensions, we take the 
same lattice as in the previous example, but prepare the 
system in a mixture of Af = and Af = 1, Af z = states. 
Recalling that for this geometry, both N s and N tz are 
conserved quantities, we find the constraint on the varia- 
tional wave function ipi X = tpi y = and ipi Z = e lipi . With 
this constraint, the dipole interaction energy is 



E, 



dip 



7 sin 2 (29) ^-^ cos(ipi — ipj) 
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(8) 



Here, the dipoles are confined to point in the z-direction, 
and therefore cannot point head-to-tail. This gives an- 
tifcrromagnctic ordering in all directions, q = (n, it, 0), 
with ipi Z = e l ( q ' Ri+< w where ipo is an arbitrary phase. 
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FIG. 1: The ordering wave vector as the lattice is tilted by 
angle a. As described in the text, for this situation a magnetic 
field is used to break the degeneracy between (jV = 1, N z = 
— 1, 0, 1) states. 

For the final example for the Mott insulating state, 
we consider a lattice in the xz-plane given by ai = 
cos(a)x + sin(a)z and a 2 = sin(a)x + cos(a)z. In addi- 
tion, we consider breaking the degeneracy of the (N = 1 , 
Af z = —1,0,1) states with an external static magnetic 



field in the z-direction which will introduce the term pro- 
portional to BL Z into our hamiltonian. Preparing the 
system in a superposition of Af = and Af = l,Af z = 1 
states, we note that because of this degeneracy breaking, 
there will be no mixing between other angular momen- 
tum states. That is, we can completely neglect the t_i ; o 
states. This will give tpi X = —iipi y = e lfi /^/2 and "02=0 
which will confine our dipoles to rotate in the xy-planc 
as: (di(t)) = doCOs(ipi — 2B e t)x + dosm(ipi — 2B e t)y. The 
dipolar energy of this system is therefore 



7sin 2 (6») ^ cos(<^ - ipj)(l 
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We use the ansatz tpi = q ■ R, + ipo to find the minimum 
of this dipolar energy for a particular lattice defined by 
the angle a, and the results are summarized in Fig. ^ 

We now consider melting the Mott insulator, and en- 
tering the superfluid (SF) state. An interesting question 
to consider is what happens to the ordering wave vec- 
tor as the Mott insulating state is melted? For instance, 
deep in the superfluid phase, we will have q — which 
is favorable for Bose-Einstein condensation while we saw 
that antifcrromagnctic ordering is typically favored in the 
Mott insulating state by dipolar interactions. One pos- 
sibility is that the wave vector interpolates smoothly be- 
tween these two extremes as the hopping J increases. An- 
other possibility is that the molecules in the s and t states 
phase-separate. We will show below that both scenarios 
are possible depending on on-site energy parameters in 
our original hamiltonian. For simplicity, we restrict our 
attention to the third example we discussed above for 
the Mott insulating state which was a two dimensional 
lattice in the xy plane prepared with a z polarized light. 
For further simplicity, we take (N s ) = (N t _) = N/2. As 
we saw before, wc can neglect populating the t x and t y 
states, and this phase has antifcrromagnctic q = (it, w, 0) 
order in the Mott insulating phase. 

Allowing for nonintcger occupation per site motivates 
the variational wave function 



i*>-n(i>$> 



(10) 



where a\ = cos(#)e* Pa ' Ri sj +sin(6*)e ip «' R< t]^ and normal- 
ization requires ^ n \a n \ 2 = 1 (compare with ©). As 
before, this wave function maximizes the dipole energy 
for a given site which is energetically favorable. We can 
now use a canonical transformation to write our origi- 
nal hamiltonian in terms of the boson operators a\ (de- 
fined above) and b\ = - sin(6»)e ip - R -s l t + cos(6»)e i P' R "4 
(a new variable resulting from the transformation), and 
drop the terms which give zero when evaluated using the 
above variational wave function (|10fl . This leads to the 
following single-site mean field hamiltonian 
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H U f = -2J J2 y cos4 ( d ) +sin 4 (6») + 2cos 2 (6»)sin 2 (6»)cos( gQ ) (a* (a) + a(a f ) - (a f )(a)) 

a— x,y 

+ Jsm 2 (20)(27i a (n a )-(n a ) 2 ) ^ 



(11) 



a—x,y 

7 . 2/ om /r, / \ , , 2 ^r cos(q-R,) 1 1 



+ 2^n»(n a - 1) + ^ (V - (7) sin 2 (20)n a (n a - 1) 
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FIG. 2: The phase diagram for V = (7 (left) and V = 
(right). For both cases, the dipolar interaction strength was 
fixed at 7 = U /5. Shown are the antiferromagnetic MI states 
with one and two bosons per site labeled Mil and MI2. SF1 
and SF2 correspond to superfluid states with partial and com- 
plete phase separation (described in text). SF3 is a superfluid 
phase with no phase separation which has an ordering wave 
vector that interpolates between the Mott insulating and deep 
superfluid regime. 



where n a = a) a and we have already performed the 
minimization over the center of mass momentum p = 
(Pt + Ps)/2. The ground state of this hamiltonian for 
fixed 9 (relative concentrations) and q = p< p s (rela- 
tive momentum) can be determined self-consistently in 
(a) and (n a ) through iteration numerically. The general 
approach will then be to minimize these ground state 
energies over q X:V <E [0, 7r] and 9 € [0, n/2]. When the 
minimum occurs for 9 7r/4, phase separation will oc- 
cur. 

The resulting phase diagrams are shown in Fig. EI The 
Mott insulating phases are antiferromagnetically aligned 
and were discussed in the previous section. SF1 cor- 
responds to partial phase separation where part of the 
lattice will have a larger concentration of s molecules 
while the other part will have a higher concentration of 
t molecules. Recall that phase separation will occur for 
when 9 7^ 7r/4 since we initially prepare the system to 
have equal populations of molecules in the s and t z states. 
The region with more s molecules will have (p s )x,y = 
and (j>t)x.y = ft- This will allow the more populated s 
species to benefit maximally from BEC which prefers zero 
wave vector while still giving q XtV = ir which is preferred 
for the dipole interaction. The similar situation holds for 
the region of the lattice with a higher concentration of 
t z molecules. SF2 corresponds to the case where the s 
and t z molecules completely phase separate. Since the 
dipole interaction is negligible for this case, we will have 



(Ps)x,y = (Pt)x,y = which will favor BEC. Finally, SF3 
corresponds to the case mentioned above where the wave 
vector q interpolates between the deep superfluid and 
Mott insulating states (0 < q XlV < ir) for which no phase 
separation occurs (9 = 7r/4). 

In conclusion, we have shown that polar molecules pre- 
pared in a mixture of two rotational states can exhibit 
long-range dipolar interactions in the absence of an exter- 
nal electric field. We have described several novel Mott 
insulating and superfluid phases that can be realized as a 
result of such an interaction. Such states can be detected 
by Bragg scattering or by time-of- flight expansion |25j . 

This work was supported by the NSF grant DMR- 
0132874 and career program, the Harvard-MIT CUA, 
and the Harvard-Smithsonian ITAMP. We thank B. 
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Note added: When this manuscript was close to com- 
pletion we became aware of a paper considering a similar 
system [l^ . 
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